Minimal models of weighted scale-free networks 
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We consider a class of simple, non-trivial models of evolving weighted scale-free networks. The 
network evolution in these models is determined by attachment of new vertices to ends of preferen- 
tially chosen weighted edges and by updating the weights of these edges. Resulting networks have 
scale-free distributions of the edge weight, of the vertex degree, and of the vertex strength. We dis- 
cuss situations where this mechanism operates. Apart of stochastic models of weighted networks, we 
introduce a wide class of deterministic, scale-free, weighted graphs with the small-world effect. We 
show also how one can easily construct an equilibrium weighted network by using a generalization 
of the configuration model. 

PACS numbers: 05.50.-|-q, 05.10.-a, 05.40.-a, 87.18. Sn 



I. INTRODUCTION 

All real-world networks are weighted. Let us explain 
this strong claim in more detail. Usual objects of in- 
terest in the science of networks are relatively simple 
nets where all edges have equal "weights" li, J., _3, 4,, 5] . 
The elements of the adjacency matrices of these net- 
works are ones and zeros. We stress that these simple 
graphs are only parodies of real networks where con- 
nections between vertices are not equal. Edge weights 
are introduced to describe this diversity. In the simplest 
case, which we discuss here, wei ght s are positive numbers 

p,7','?,'?,'To','iiEiiiiaiiiMiiiiaiiiii,i2ai2a 

(see Fig.0. So, the elements of the resulting adjacency 
matrices are zeros and positive numbers, Wij. However, 
in principle, in more complex situation, edges may be de- 
scribed by a set of variables or operators. Note that an 
extra set of variables may be introduced to describe an 
individual properties of vertices (so called hidden vari- 
ables 23, 24, 25, 21 mm ). 

Edge weights allow one to better describe reality. For 
example, in the simplest, unweighted version of a network 
of social contacts, each edge connects a pair of individuals 
who had one or two or more contacts. In a far more infor- 
mative weighted network of social contacts, the weight of 
an edge shows the frequency of recent contacts between 
the corresponding pair of persons l29l l30| . Another ex- 
ample is a network of coauthorships, where edge wei ght s 
show the number of joint papers of two coauthors 
(This network may be treated as a one-mode projection 
of a bipartite graph.) Note that networks with multi- 
ple connections can be considered as weighted ones with 
integer edge weights. 

Weighted networks provide an informative, usable pic- 
ture of reality. Our first example is collaboration net- 
works. Collaborations are adequately represented by bi- 



partite collaboration graphs, where one kind of vertices 
is collaborators and the other one is acts of collaboration 
[see Fig. |3 graph A] . Usually, empirical researchers con- 
sider one-mode projections of these graphs, where a pair 
of vertices-collaborators is interconnected by an edge if 
there has been at least one act of collaboration between 
them [see Fig.|21 graph B]. This net is more simple than 
the original bipartite graph, but an essential information 
is lost. In particular, one cannot recover the bipartite 
original by using this projection. The weighted one-mode 
projection shown in Fig. 13 graph C, is more informative 
than unweighted projection B but, nonetheless, is less in- 
formative than original bipartite graph A. One still can- 
not recover bipartite graph A if weighted projection C is 
known. 

The second example is energy landscape networks, 
where each vertex is a local minimum of the potential en- 
ergy landscape, and an edge connects two minima with 
a saddle point between them |^ 1^3 • An unweighted 
version of this network reflects only a general structure 
of connections in the configuration space of a system. In 
contrast, directed weighted graphs, where edge weights 
indicate transition rates provide basic necessary informa- 
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FIG. 1: (a) A typical one-partite, weighted, 
graph, (b) A bipartite weighted graph. 
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FIG. 2: One-mode projections of a bipartite collaboration 
graph — graph A. The circles and squares show the collabora- 
tors and the acts of collaboration respectively. Graph B is the 
unweighted one-mode projection of graph A. It is impossible 
to recover bipartite graph A by knowing unweighted graph B. 
Graph C, which is the weighted one-mode projection of graph 
A, is more informative than graph B but, nevertheless, one 
still cannot recover graph A if the weighted projection C is 
known. Edge weights show the number of collaboration acts 
between the corresponding vertices. 

tion about the relaxation and kinetics of a system. (One 
may also ascribe the energies of the potential minima to 
the vertices, and the energies of the saddle points to the 
edges.) 

The third example is various chemical reaction net- 
works (e.g., networks of metabohc reactions [131 )■ In 
a one-partite version of these networks, weights show the 
chemical reaction weights. Networks of corporate own- 
erships are also directed weighted graphs (edge weights 
indicate the fractions of company's shares in hands of 
other companies) . The links of spacial networks [^3 has 
a natural characteristic — their length. So, in this case an 
edge weight may be introduced as some function of the 
length of the corresponding edge. 

The simplest local characteristics of weighted networks 
are (i) the weight of an edge, u;^, (ii) the degree of a 
vertex, and (iii) "the strength" of a vertex, 

S, EE W,j , (1) 

Here the sum is over the nearest neighbors of the ver- 
tex i. (One can introduce directed weighted networks, 
where vertices have "in-strength" and "out-strength".) 
The distributions of these local characteristics in ran- 
dom networks are a weight distribution Q{w), a degree 
distribution P{k), and a strength distribution R{s). 

A number of models, which provide complex, in par- 
ticular, scale-free weighted networks, were proposed (see, 
e.g., Refs. IE mH 111 111 El 0)- Obviously, one may 
easily construct a complex weighted network indepen- 
dently ascribing degrees to vertices and weights to edges 
(see the next section). For this, (i) build an unweighted 
network by any of existing algorithms and afterwords (ii) 



arrange weights of edges. These are trivial constructions. 
A much more interesting scale-free network model, where 
the evolutions of degrees and weights are coupled, was 
proposed in Refs. |l3|ll','T^. 

In this network, in addition to a standard preferential 
attachment rule HHs^ (see also H^), a specific redis- 
tribution of edge weights was introduced. In more detail, 
(i) each new vertex is attached to a preferentially se- 
lected existing vertex (the probability of the attachment 
is proportional to the vertex strength), and (ii) weights 
of connections of the latter vertex are updated in a spe- 
cific way depending on the strength of the vertex. That 
is, in this model, a preferential attachment to a vertex 
induces changes in weights of its edges. This evolution 
results in (a) a power-law weight distribution, (b) a linear 
dependence of the strength of a vertex on its degree, and 
(c) power-law strength and degree distributions. Note 
that the specific form of rule (ii) is necessary to obtain 
scale-free networks. 

In the present paper we show that quite similar results 
[(a), (b), and (c)] can be obtained by using a more simple 
construction. Actually we use the same idea as in our 
network evolving due to attachment to the ends (or to 
an end) of a randomly chosen edge Is^. In the present 
network, 

(i) the weight of a preferentially chosen edge is in- 
creased by a constant, and 

(ii) its end receives a new connection of a unit weight. 

One can see that this model is so simple, that it may be 
considered as a minimal one. 

The paper is organized as follows. In Sec. |n] we, for 
comparison, demonstrate a "trivial" construction of an 
equilibrium weighted network, which is a direct general- 
ization of the configuration model. In Sec. IIIII we define 
our stochastic model and present its solution. In Sec. IIVI 
we explain how one can choose edges at random or ac- 
counting for their weights. In Sec. we describe simple 
generalizations of the model. In Sec. I VII we introduce de- 
terministic scale-free weighted networks. In Sec. IVIll we 
discuss possible applications of our models. 

II. EQUILIBRIUM WEIGHTED NETWORKS 

As a starting point, we consider a "trivial" equilib- 
rium, weighted network — the simplest generalization of 
the configuration model of a random network. The stan- 
dard configuration model is, in simple terms, a maximally 
random graph with a given degree distribution, -P(fc), see 
Ref. where the model was introduced in this form, 
and also Refs. [4lll4^l4^ with very similar constructions. 
In addition, we assume that the weights of the edges in 
the random graph ensemble of the configuration model 
are random, independent, and distributed according to a 
given weight distribution Q{w). The question is: what is 
the resulting distribution of the vertex strength, i?(s), in 
this network? 
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For brevity, in this section, we consider only networks 
with integer edge weights. However, the final result for 
the asymptotics of distributions is also valid without this 
assumption. 

Actually, we must calculate the distribution of the 
sum of independent equally distributed random vari- 
ables. The resulting relation is especially simple in a 
Z-transformed form. Let us introduce the Z-transforms 
(generating functions) of the distributions P{k), Q{w), 
and i?(s): 

$(.)^^.'=P(fc) 
fc 

I](z)^^z-Q(u;) 

W 

vl/(z)^^z«i?(s). (2) 

s 

Then, accounting for the definition of the strength 
and for the mutual independence of edge weights (and 
vertex degrees), we immediately have 

vI'(z)===^P(fc)S'=(z) = ci>(E(z)). (3) 

k 

Here we used a standard relation for the Z-transform 
of the distribution of the sum of independent random 
variables (see, e.g., Ref. 13). 

The asymptotic behavior of the original distribution at 
large values of a random variable is related to the ana- 
lytical structure of the corresponding Z-transformation 
near z = 1. For example, for a power-law original, we 
have the following correspondence: 

P{k > 1) ~ k^^" 

^'(znearl) = 1 + ^ a,(l - z)* + 6(1 - z)'"'-^ (4) 

i>0 

On the other hand, the ^-transformation of a rapidly 
decreasing distribution (i.e., with the finite moments), is 
analytical at z = 1. So, accounting for correspondence 
(0J) , relation ^ leads to the following conclusions for this 
network: 

(i) if both the distributions P{k) and Q{w) are rapidly 
decreasing, then the strength distribution R{s) is 
rapidly decreasing too; 

(ii) if P{k) is rapidly decreasing, and Q{w) is, e.g., 
a power law, Q{w) ~ w~'^"' , then the resulting 
strength distribution is also power-law, with the 
same exponent, R{s) ^ s~'^"' , i.e., 7^ =7™; 

(iii) if P{k) is, e.g., scale-free, P{k) ^ fc~'"=, and Q{w) 
rapidly decreasing (i.e., S(z) analytic at z = 1), 
then R{s) s"'^^ i.e., 7^ = 7^; 




FIG. 3: Schematic view of the network growth. At each time 
step the weight of a preferentially chosen edge is increased by 
a constant A, and a new vertex is attached to the ends of this 
edge. 

(iv) if both P{k) and Q{w) are scale-free, then the re- 
sulting distribution R{s) is also scale-free with ex- 
ponent 7s = min(7fe,7„). 

These are important, although particular, cases. How- 
ever, relation (j^J allows one to easily obtain, after 
an inverse Z-transform, a strength distribution for the 
weighted graph with arbitrary degree and weight distri- 
butions. 

Note that, also, one may consider a more "serious", 
"non-trivial" generalization of the configuration model. 
These are labeled random graphs with a given sequence 
of "generalized degrees". The generalized degree of a 
vertex is a complete set of numbers, where each number 
shows how many edges of a given weight are attached to 
this vertex. This is quite similar to the constructions of 
graphs with hidden color of Refs. [isf. 

III. ATTACHMENT TO WEIGHTED EDGES 

Let us introduce the model of a growing scale-free net- 
work in more precise terms than in the Introduction. We 
first discuss the simplest case. 

We assume that the growth starts from an arbitrary 
configuration of vertices and edges, e.g., from a single 
edge of weight 1. At each successive time step, 

(i) choose an edge with probability proportional to its 
weight and increase this weight by a constant A > 
0, 

(ii) attach a new vertex to both the ends of this edge 
by edges of weight 1 

(see Fig.EJ. 

In the particular case of A = 0, this network is reduced 
to our model of Ref. [s^. Rule (i) may be interpreted 
in the following way. Suppose that A and the weights 
of the edges are integer numbers. In this case, an edge 
of weight w may be treated as a w-multiple edge. If 
we choose each "elementary" edge in this network with 
multiple edges at random, with equal probability, we just 
arrive at the proportional preferential choice of edges in 
the corresponding weighted network. 

Rule (ii) may be easily modified: attach a new vertex 
to one of the ends of the selected edge and not to the 
both of them — see below. In the present form, rule (ii) 
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results in the following rigid coupling of the degree ki of 
a vertex and its strength sf. 



fc,(l 



2A. 



(5) 



Indeed, in this model, each attachment to a vertex in- 
creases its degree by 1 and its strength by 1 + A (1 — due 
to the link of weight 1 to a new vertex and A due to 
the modification of the selected edge). Taking into ac- 
count that new vertices have degree and strength equal 
to 2 results in relation (0). So the degree and strength 
distributions have quite similar asymptotic behaviors. 

We will show below that the weight distribution, the 
strength distribution, and the degree distribution of this 
network have power-law asymptotics: 



QH-w-"^-, i?(s)-s-'^% P{k)^k-'"' (6) 
with exponents 



In 
Is 



2 + A 2 
A A 
, 2 + A ^ 

^^ = ^ + TTa^2- 



1 



1 + A 



(7) 

(8) 



Formula (TJ is valid at A > 0, and formula © is valid at 
A > 0. 

Note that formulas 0, © give 7s = 7fe < 7u.. Recall 
the relation 7^ = min(7fe, jw), which we have obtained for 
scale- free equilibrium networks (the "configuration model 
of a weighted network" ) . So, formulas Q, © satisfy that 
relation. One can easily check that the same is valid for 
the other scale-free weighted networks in this paper [see 
formulas (jHI and (HHJl, (CHI) and 1(201), lEU and 

Note that the evolution of edge weights in this model is 
quite similar to that in the Simon model (4Sj . One must 
take into account that in our network the total weight 
of edges growth proportionally to t: W(t) = t{2 + A), 
which is asymptotic expression at large t. This gives 
the following evolution equation for the mean number 
N{'w, t) of edges of weight w at time t: 



N{w,t + l) = N{w,t) +2dn,,i 



t{2 + A)' 



t{2 + A) 



N{w,t). (9) 



(For similar equations for degree distributions in growing 
scale-free networks, see Refs. ^3 and At each time 

step two new edges of weight 1 emerge. This produces 
the term 25^,^1 on the right-hand side of the equation. 
The third term describes increase in the number of edges 
of weight w due to the modification of edges of weight 
w — A. The last term on the right-hand side describes 
the reduction of the number of edges of weight w due to 
their modification. 



We pass to the continuum t limit and assume a 
stationary form of the weight distribution, N{w,t) = 
2tQ{w,t) — > 2tQ{w). Then, passing to the continuum 
limit of weight allows us to obtain a power-law form of 
the weight distribution Q(w) ~ w~'^"' with the 7^ expo- 
nent given by formula l(7|l. 

The exact solution of the stationary limit of Eq. Hl()|l , 



Q(w) cx B{w/A,j„ 



(10) 



is quite similar to the solution of the Simon model (we 
have omitted a normalization factor on the right-hand 
side of this formula). Here B{ , ) is the beta-function. 
As is natural, at large w, this expression has a power-law 
form with the same exponent given by relation |7J. 

The total strength of the vertices in the network is 
S{t) = 2W{t) = 2t{2 + A). [The total degree of the 
network evolves in the following way: Kit) = At.] The 
attachment of a new vertex to the end vertices of an edge 
increases their strengths by 1-1- A. Furthermore, the rules 
of the model actually result in preferential (proportional) 
attachment to vertices of higher strength. So, the evolu- 
tion equation for the mean number N{s, t) of vertices of 
strength s at time t looks as follows: 



iV(s,t + 1) = N{s,t) + 8, 



s,2 



^(1 + A) 



N{s~(l + A),t) 



-N{s,t). (11) 



2^(2 -f A) ' V ' 2t{2 + Ay 

The factors 2 of the third and the fourth terms on the 
right-hand side of this equation are due to the attachment 
of each new vertex to the two ends of a selected edge. The 
resulting stationary strength distribution R{s), N{s, t) = 
tR{s, t) R{s), is of the form 



i?(s) cxB(s/(1-HA),7, 



(12) 



with 7s given by formula ^ . Consequently, due to linear 
relation © between strength and degree, the resulting 
degree distribution is 

P(fc)^B(fc,7fc) (13) 

with 7fc given by formula © . One can see that the result- 
ing expressions of exponents, Eqs. lO andj(7|^, are similar 
to those obtained for the model of Refs. ,1(1 . .ill fl^ . The 
ranges of variation with the parameter of a problem — 
(2, 00) for 7u, and (3, 00) for 7s and jk — are the same in 
both these networks. 

By construction, our networks have numerous loops 
of length three, which indicates high clustering (see 
Ref. [33 )■ However, one should be careful with this 
claim. There are two distinct integrated clustering char- 
acteristics. One can show that the average clustering (the 
average clustering coefficient of a vertex) in this case is 
finite in the infinite network limit. On the other hand, 
the clustering coefficient (in simple terms, the density of 
triangles in the network) approaches zero in the infinite 
network. 
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IV. HOW TO CHOOSE AN EDGE AT RANDOM 

A meticulous reader may say: "Your model, as it was 
formulated, is indeed looks more simple than that of 
Refs. OdUilLJ. However, your model is based on the 
selection of connections, which may be, computationally, 
more complex task than the selection of vertices. This 
may depreciate this model." 

To convince this reader, we indicate two easy (and 
quick) ways to select connections, (i) Make a list of 
edges (this list is sufficiently short in sparse networks) 
and choose from it. (u) Choose edges, starting the proce- 
dure from random vertices. Below we describe the second 
way in more detail. 

Let us first demonstrate how one can select at ran- 
dom edges in a network. The simple procedure looks as 
follows: 

(i) choose a random vertex and then 

(ii) subsequently choose each of its edges with some 
sufficiently small probability p, 

then, repeat. Here, the probability p must be small: 
p < 1/fcniax, or at least p -C 1/fc, where fcmax is the maxi- 
mal degree of a vertex in the network, and k is the mean 
degree. This condition is necessary to subsequently se- 
lect single edges but not bunches of edges — "hedgehogs". 
(Note that, unlike the above procedure, if we simply se- 
lect a random edge of a random vertex, we will not choose 
a random edge.) 

Similarly, in a weighted network, one can choose an 
edge with probability proportional to its weight. For this, 

(i) Choose a random vertex and then 

(ii) subsequently choose each of its (weighted) edges 
with some sufficiently small probability propor- 
tional to its weight. 

then, repeat. Rule (ii) means that, after a vertex (say 
vertex i) is selected, each of its weighted edges is selected 
with probability pwy /sj. Here, Si is the strength of this 
vertex [see definition J^l] and Wij is the strength of an 
edge attached to this vertex, p is an arbitrary parame- 
ter, so small that there is a small chance to subsequently 
select two or more edges attached to one vertex (see the 
above procedure). 



SIMPLE GENERALIZATIONS 




FIG. 4: The scheme of the growth process for the modified 
network. At each time step the weight of a preferentially 
chosen edge is increased by a constant A, and a new vertex 
is attached to any one of the ends of this edge. 



(ii) attach a new vertex to one the ends of this edge by 
edges of weight 1 

(see Fig.lU. 

The total weight, the total degree, and the total 
strength of the vertices of this network grow in the follow- 
ing way: W{t) ^ t(l + A), K{t) ^ 2t, S{t) ^ 2t{l + A), 
respectively. We will show that the network has power- 
law distributions (j^ with exponents 



1^ = 1 + = 2 + -r , 



1 + A 



1-I-2A 



(14) 
(15) 



Instead of rigid coupling between the degree and 
strength of an individual vertex in the network of Sec. lHIl 
now we have only the asymptotic relation between the 
mean values of the strength and degree of an individual 
vertex: 



{l + 2A)ki 



(16) 



This relation is valid for highly connected vertices (with 
large degrees). It shows that 7s = 7fe in this network. 

Asymptotic equality H16(l follows from the following 
simple considerations. At each modification of vertex i, 
its degree and strength are modified in the following way: 
(i) with probability 1 /2, ki —>■ ki + 1 and + A + 1, 

and (ii) with probability 1/2, ki — > ki and Si Si + A. 
This directly leads to relation (fTH| . 

One can see that for this network, the equation for the 
mean number of edges of weight w at time t is 



As we promised, let us consider the variation of the 
model with attachment to one end vertex of a selected 
edge. At each successive time step, w — A 



N{w,t+1) = N{w,t) +dyj,l 



■N{w - A,t) - —N{w,t) (17) 



t{l + A) ' ' ' t{l + A) 

3 with Eq. lO]. This immediate! 
0, mula H14|) for the 7^ exponent of the weight distribution. 



(i) choose an edge with probability proportional to its 

weight and increase this weight by a constant A > [compare with Eq. I^] . This immediately leads to for- 
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For the mean number of vertices of strength s in this 
network, we have the foUowing equation: 



W[+ A 



N{s,t+1) = N{s,t) + Ss 



»-<' + ^'«(,-(l + A),,)- 



+ 



2t(l + A) 

s-A 
2i(l + A) 



' N{s,t) 



2t{l + A) 



N{s - A,t) 



2t{i + a; 



(18) 




W2+ A 



w,+ A 



[compare with Eq. . This, together with relation 
leads to formula 1151) for the 7^ exponent of the strength 
distribution and the "fk exponent of the degree distri- 
bution. One can see that the only difference between 
formulas 0, ^ and (fnji . l(T5|l is that the parameter A 
in the results of Sec. Illll is now substituted by 2S. 

In the same way as above, one can consider a combina- 
tion of the models that we have already discussed. For ex- 
ample, the following combination of the evolution chan- 
nels may be introduced. At each time step: (i) Choose 
preferentially an edge and increase its weight, (ii) (a) 
With probability p2 , attach a new vertex to both the ends 
of this edges; (b) with probability pi , attach a new vertex 
to one of the end vertices of this edge; (c) with probabil- 
ity pq, do not add a new vertex. Here, Pq + Pi + P2 — 1- 
One can also add the processes of linking of existing ver- 
tices, rewiring of connections, and introduce variation of 
A. 

Up to now we discussed the proportional preferential 
choice of edges, that is, edges were chosen with prob- 
ability proportional to their weights. In general, one 
may choose edges with probability proportional to some 
function j(w) of their weight (preference function). Fur- 
thermore, one may consider an inhomogeneous situa- 
tion, where a preference function depends on an edge 
fij{w). Here the preference functions fij{w) (and their 
parameters) may be randomly distributed similarly to 
what was considered for more traditional preferential at- 
tachments to vertices in inhomogeneous networks (see 
Refs. HOl)- A also may be made random variable. 
One may also, at each time step, attach a new vertex to 
ends of several edges but not to a single edge as above 
(see Fig.|S}. By using these generalizations, we can easily 
obtain various complex weighted network architectures 
with fat-tailed and rapidly decreasing distributions, with 
condensation, gelation, etc., previously studied in un- 
weighted networks (42j \iSk 1^ 1^ ■ 

Let us suppose, for example, that (i) each new vertex 
becomes attached to all the ends of m > 1 preferentially 
selected edges, (ii) the preference function is a linear func- 
tion: f{w) = w + a where a > — 1, (iii) weights of the 
selected edges are increased by a constant A. In this 
case, formulas O and ijSJl take the forms: 



FIG. 5: The growth process with attachment of a new vertex 
to several weighted edges. 



w 




FIG. 6: A generalization of the network growth process 
shown in Fig. |21 The weights of new edges, Ai and A2, and 
the addition A3 may be random. 



7«> = 2 + 2 



1 + a 



A ' 



7s = 7*; = 2 + 



1 + a 



1 + A + a 



(19) 
(20) 



The derivation of these expressions is very similar to that 
in Sec. IIIII [one must take into account relation If 
a —1, then all the three exponents ^w,Js,Jk ^ 2. 
Note that if a ^ 00, then 7s,7fc —>■ 3. One can easily 
understand this limit. The infinite a actually means the 
absence of preference, and edges are chosen at random, 
without accounting for their weights, exactly as in our 
network of Rcf. 39]. So, we arrive at the same value of 
exponents as in that model. (Recall that 7^ = 3 in the 
Barabasi- Albert model ^3^.) 

Another natural generalization of the process in Fig. |3| 
is shown in Fig. |B| In this process, new edges also can 
have different weights. Furthermore, all three numbers, 
Ai, A2, and A3, may be random. 



VI. DETERMINISTIC WEIGHTED 
NETWORKS (PSEUDOFRACTALS) 

Compact growing networks with the small-world ef- 
fect m ay b e produced in a deterministic way 
HI El El El [58, 59, 60]. These graphs have a discrete 
spectrum of degrees. These spectra may have a variety of 
shapes: they may be fat-tailed, in particular, scale-free. 



they may be rapidly decreasing, e.g., exponential. The 
key and necessary feature of these deterministic graphs is 
the small- world effect: their mean intervertex distances 
grow slower than any (positive) power of the numbers of 
vertices. So, while these graphs visually look very similar 
to fractals, they are only parodies of fractals. The dif- 
ference is crucial. Fractals have a finite dimension (i.e., 
their intervertex distances grow as a power of the num- 
bers of their vertices). In contrast, the networks, which 
we discuss, have infinite dimension (i.e., their intervertex 
distances grow slower than any power of the numbers of 
vertices, e.g., logarithmically). This is why we call deter- 
ministic graphs of this kind pseudofractals [3. l52l| . 

One should stress that a scale-free architecture (a 
power-law degree spectrum) is not not a definitive, nec- 
essary feature of these graphs. Many fractals have scale- 
free spectra of degrees (see discussion in Ref. Q). 

How may these deterministic compact networks with 
a self-similar structure be constructed? There are two 
ways to obtain a compact deterministic network: 

(i) The first was realized in Ref. [sj. At each step of 
the evolution, generate a number of the copies of a 
graph and connect them together, e.g., by adding 
new edges, which decreases intervertex distances. 

(ii) One can use the following approach [2, E^- At 
each step, in a regular manner, transform given ba- 
sic subgraphs of the network into other, larger but 
compact, configurations consisting of the same ba- 
sic clusters. These basic clusters may be vertices, 
edges, triangles, squares, etc. The words "trans- 
form in a regular manner", in the simplest, very 
particular case, mean "transform all the given ba- 
sic subgraphs of the graph". The words "larger 
but compact configurations" are illustrated by an 
example below. 

(Note that in some specific situations, these approaches 
are equivalent.) 

Let us demonstrate how to construct a deterministic, 
scale-free, weighted, compact network by using the sec- 
ond approach. We will exploit the equivalent represen- 
tation of integer-weighted edges as multiple connections 
(see Fig. 13). 

As an example, we use the transformation that is a 
deterministic variation of the transformation shown in 
Fig. |31 with an integer A. This transformation is shown in 
Fig. |S1 The edge of integer weight w is transformed into 
the edge of weight w{l + A) with w triangles of 1-weight 
edges attached. The growth starts from the triangle of 
edges of weight 1. At each successive step, each edge 
of the graph is transformed in the way shown in Fig. |S1 



w=4 



(a) 




(b) 




(c) 




FIG. 8: (a) An example of transformation of "elementary" 
edges. A = 1 (compare with Fig. I^J. (b) An equivalent rep- 
resentation of the transformation of an edge of weight 1. (c) 
The resulting transformation of an edge of weight 2 following 
from (a) or (b). 



The result is a weighted pseudofractal network, shown in 
Fig- m in particular case of A = I. If A = 0, we get the 
deterministic graph of Refs. |3,I13- 

We will show below that the weight, strength, and de- 
gree distributions of this graph are power-law with expo- 
nents 



ln(3 + A) 
ln(l + A) 



Is = 7k 



1 



ln(3 + A) 
ln(2 + A) ' 



(21) 



(22) 



respectively. Note that 7u),7s,7fe ^2 as A ^ oo. 

We emphasize that this is only a particular example. 
The readers can easily consider numerous variations: use 
other transformations of this kind or their combinations, 
transform different clusters, transform not each of given 
elements of a graph but only some of them, selected in 
a regular way, consider networks embedded in a finite 
dimension space, consider trees, and so on. These varia- 
tions produce a wide range of network architectures and 
a wide range of the exponents values in the range (2, oo). 

Let us consider the growth of the graph. One can see 
that at each step of the evolution, a weight of each in- 
dividual edge of the graph, a strength of each individual 
vertex and its degree transform in the following way: 



w —f w' — {1 + A)w . 
s ^ s' = (2 + A)s , 



FIG. 7: Multiple edges of unweighted networks may be tre- 
ated as single weighted edges. 



(23) 

respectively (the indices of edges and vertices are not 
shown) . 



8 



n=0 



n=l 



n=2 




FIG. 9: An example of a deterministic, scale-free, weighted 
graph with the small- world effect. The graph is obtained by 
successive application of the transformation from Fig.|H|(A — 
1) to all the edges of the graph. The evolution starts from the 
triangle of edges of weight 1 (n = 0) . Only the first two steps 
are shown. The bare edges denote the edges of weight 1. 

We introduce the following notations for basic num- 
bers: Nn, Ln, and Wn, which are the total numbers of 
vertices and edges and the total weight of the edges in 
this deterministic graph in an nth generation, respec- 
tively. The total degree is Kn = 2_L„. The total strength 
of vertices is Sn = 2Wn- 

In the initial state, Nq = Lq = Wq = 3. Then, one can 
obtain 

Wn = 3(3 + A)" . (24) 

One can see that 

L„+i = L„ + 2Wn . (25) 

This gives 



in = ^^[2(3 + A)" + A]. (26) 

To find the weight, strength, and degree spectra we 
take into account the following circumstances. New edges 
have weight 1, new vertices have degree and strength 
equal to 2. All the edges that emerge simultaneously have 
the same weight. All the vertices that emerge simultane- 
ously have the same strength and the same degree. The 
weight Wm of edges emerged m < n generation ago and 
the strength Sm and the degree km of vertices emerged 
at that moment are 

Wm = (1 + A)™-i , 
= 2(2 -f A)'"-i , 

km = Y^ii^ + ^'r''+^]- (27) 

The nth generation graph contains Wn-\ — 3(3 -I- A)""-'^ 
new vertices and 2Wn-i = 6(3 -I- A)"~^ new edges. 
Consequently, the nth generation graph contains 6(3-1- 
A)"-™-i edges of weight Wm and 3(3 + A)"-"-i ver- 
tices of strength Sm and deg rcG kfji. Here Wm-) 

and 

km are given by formula H27I) . 

As a result, we obtain a power-law behavior of the spec- 
trum of the weight: the number of the edges of weight Wm 
is N{'Wm) oc w^'"*''^^^-'/'"'^^^''. This is a discrete spec- 
trum with gaps between discrete weights growing with w. 
To present our result in a form which can be compared 
with results for stochastic models, where distributions 
are continuous and spectrum gaps are absent, we pass 
to the cumulative weight distribution: X]iu,>iu N^wi) oc 
Wm'"^'^^'^''^ in(i+A) ^ r^Yiis distribution already has no gaps 
like the corresponding cumulative distribution of the edge 
weight of stochastic networks, Ncumiw) oc w"*^^™"^'. So 
we arrive at formula (|21|l for the jw exponent of the edge- 
weight distribution. 

Similarly, we obtain the power strength spectrum: 
the number of vertices of strength Sm is N(sm) oc 
Sm'"''^^'^''^'"''^^'^''. The cumulative distribution of the 
vertex strength is J2si> ^i^i) oc and 
the corresponding cumulative continuum distribution of 
the vertex strength in stochastic models is Ncuni{s) oc 
g-(7s-i)_ One can see that the degree spectrum of the 
deterministic graph is quite similar to the strength spec- 
trum. So, we finally obtain formula (|22|l for the expo- 
nents of the strength and degree distributions. 

Other structural characteristics of this graph also can 
be easily found. We present here only an expression for 
the standard local degree-dependent clustering, which is 
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defined as 



Cik): 



mean # of triangles attached to a vertex of deg. k 



k{k -l)/2 



(28) 



By construction, a vertex of degree k in this graph has 
k — 1 triangles attached, so we readily find 



C{k) = 2/k 



(29) 



Note that this expression is independent of A, and we 
have the same result as for the network of Ref. la, |53 • 



VII. 



THE RATIONALE BEHIND THESE 
MODELS 



We discussed networks where edges of high weight at- 
tract new connections. Does this occur in the real world? 
Here we present only one illustration. 

Let us consider the evolution of a weighted network of 
scientific coauthorships, taken in a one-mode represen- 
tation. So, the vertices are authors, and the weighted 
edges are (pair-wise) coauthorships. Intensive pair-wise 
collaborations have a greater chance to attract new col- 
laborators than occasional connections. Speaking in sim- 
ple terms, if two coauthors have only a single paper, e.g., 
written 20 years ago, there is a small chance that some- 
body suddenly decide to write a paper with them. In 
contrast, fruitful joint efforts lead to new coauthorships 
much more frequently. One may say, it is the papers — 
scientific results, but not their authors that attract new 
coauthors. This is precisely the mechanism that is dis- 
cussed in the present paper. 



VIII. DISCUSSION AND SUMMARY 

One point should be stressed. We have shown that 
our approach leads to results which are close to those 



obtained in Refs. 0, 0, 0|. So that, these two ap- 
proaches are related. In the model of Refs. [lOl ITllll^ . 

"strong" vertices attract new connections and afterwards 
the weights of the edges of these vertices are specifically 
modified. In contrast, in our case, links of high weight 
increase their weights and attract new connections. That 
is, in one approach, the attachment is to ("strong") ver- 
tices and in the other approach, the attachment is to 
( "heavy" ) edges. This is a principal difference which is re- 
lated to distinct real situations. We have shown that the 
introduced mechanism operates in real-world networks. 

We have demonstrated that our approach allows one to 
formulate minimal, non-trivial, solvable models of evolv- 
ing scale-free, weighted networks. We have found a num- 
ber of basic structural characteristics of these networks. 
We indicated some possible generalizations. In particu- 
lar, we have introduced a wide class of weighted deter- 
ministic graphs of a pseudofractal type and have con- 
sidered the simplest scale-free, weighted, deterministic 
graph in detail. 

In summary, we have developed an effective, simple 
approach to evolving weighted networks. This approach 
allows one to extend numerous results for well studied un- 
weighted networks to a wide class of weighted networks. 

Note added. After this paper had been prepared, works 
|i61. 62. 63] , where the linking process also is determined 
by the edge weights, have appeared in the cond-mat elec- 
tronic archive. 
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